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Investigating QCD Vacuum on the lattice 
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Investigations on the structure of QCD vacuum from first principles can be done on the lattice. The mechanism 
of confinement is an example: results from lattice on it are reviewed. 



1. Introduction 

Euclidean Feynman path integral uniquely 
identifies the ground state of a field theory (vac- 
cum)§. 

The Feynman functional integral is defined as 
the limit of ordinary integrals defined on discrete 
set of points in a four dimensional box, when the 
number of points is sent to infinity filling the box 
densely; the size of the box is then sent to infinity 
to cover the whole space time. Lattice formula- 
tion is an approssimant in this sequence. 

QCD has an UV fixed point at g = (asymp- 
totic freedom): as (7 — > the physical length scale 
^ increases in units of the lattice spacing a 

7^=0 (1) 

The density of lattice points in physical units goes 
large. On the other hand a mass gap exists in the 
theory, which makes the infinite volume limit well 
defined as a thermodynamical limit. 

Lattice is a good approximant of QCD if 



(2) 



with L the size of the lattice. 

The above argument is a strong indication that 
QCD most probably exists as a field theory in the 
constructive sense, and can be defined as a limit 
of lattice formulation as go —> 0, and L is such 
that the inequality (||) is satisfied. 

The textbook quantization of QCD is based on 
perturbation theory, and the ground state is Fock 
vacuum of quarks and gluons. 

* Partially supported by MURST and by EC, FMRX- 
CT97-0122 



Phenonienological evidence exists g that the 
QCD vacuum is not the perturbative vacuum. 
The instability of Fock's vacuum is possibly the 
origin of the non Borel summability of the per- 
turbative expansion!^. 

Perturbation theory apparently works at small 
distances but is not well defined and is unable to 
describe large distance physics. 

The most clear evidence that Fock vacuum is 
not a good approximation to QCD ground state, 
is that its elementary excitations, quarks and 
gluons, have never been observed as free parti- 
cles. This phenomenon is known as confinement 
of colour and is one of the most intriguing prop- 
erties of QCD. 

The study of the mechanism of confinement is 
an important chapter of the investigation of the 
structure of QCD vacuum on the lattice, and will 
be the object of this talk. In sect. 2 I will re- 
view the experimental evidence for confinement. 
In sect. 3 I will discuss the phenomenology of the 
deconfinement transition as observed in numeri- 
cal simulations on the lattice. This will naturally 
lead us to the idea of duality, which will be the 
object of sect. 4. Progress in understanding con- 
finement as dual superconductivity of the vacuum 
will be reviewed in sect. 5 

2. Confinement: experimental evidence. 

Confinement is defined as the absence of col- 
ored particles in asymptotic states. 

The existing experimental evidence for confine- 
ment is based on the negative result of searches of 
fractionally charged particles (quarks) in particle 
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reactions and in nature. 

The cross section aq for the inclusive produc- 
tion of quark or antiquark in the process 

p + p^ q{q) + X 

at cm. energies ~ 100 GeV has an upper hmitQ 

aq < 10"*° cm^ 

to be compared with the total cross section ctt at 
the same energy ax ~ 10~^^ cm^. 

The ratio Rq ~ (Jq/ar has then the upper 
bound 

R < IQ-^^ 

while the expected value in the absence of con- 
finement is a sizable fraction of unity. 

The negative result of the search of fractionally 
charged particles in ordinary matter by Millikan- 
like experiments gives an upper limit for the ratio 
R of the quarks abundance Uq to nucleons abun- 
dance rir, 



ii-ob = — < 10 
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corresponding to the analysis of ~ 1 gr of matter. 

In the absence of confinement the expectation 
for R, RscM, in the standard cosmological model 
is RscM - 10^12 i- Again 



R/R. 



SCM 



< m 
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A ratio smaller than 10"^''', if different from zero, 
would be too small to have a natural explanation 
in any theory. The most natural interpretation 
is then that those ratios are strictly zero, or that 
confinement is an absolute property of the vac- 
uum based on a symmetry]^. 

The question is: what symmetry of QCD vac- 
uum prevents quarks to exist as free particles? 

As for gluons, they have no such characteristic 
signature as a fractional charge, and their iden- 
tification is not clearly feasible. No experimental 
data exist on gluon confinement. We shall de- 
fine anyhow confinement as absence of any col- 
ored particle as a free particle. 



3. Deconfinement Phase transition on the 
lattice. 

QCD at finite temperature can be studied on 
the lattice. The partition function Z{T) of a field 
with action S[^] is equal to the Feynman inte- 
gral in Euclidean space, with the time direction 
extending from to 1/T and periodic bound- 
ary conditions in time for bosons, antiperiodic for 
fermions 



Z{T) 



n[r'$]e 



-/3 J d^x J^^'^dxoSE [<!>{x,xo)] 



(3) 



In lattice QCD this corresponds to having a lat- 
tice of size -/V| X Nt, with Ns > A'^r, and the 
temperature is given by the inverse of the tempo- 
ral extension Nxa {a the lattice spacing). 

The value of a in physical units depends on the 
coupling constant (/3 — 2N/ g^) via renormaliza- 
tion group 

a{(3) = -5- exp(-6o/3) 



A, 



with 

&o > 
Hence 



(asymptotic freedom) 



T 



1 



A, 



exp(6o/3) 



(4) 



As a consequence of asymptotic freedom low 
temperature (confinement) corresponds to strong 
coupling (large g) or to disorder in the language 
of statistical mechanics, high temperature corre- 
sponds to order or to weak couling. 

If confinement is due to a symmetry, it has to 
be a symmetry of the disordered phase. This nat- 
urally leads to the idea of duality, which will be 
the object of sect. 3. 

The deconfinement transition is detected on the 
lattice in pure gauge theory, by looking at the 
correlator 



V{x - v) ^ {L{x)L(y)) 



(5) 



where L{x), the Polyakov line, is the trace of the 
parallel transport along the time axis across the 
lattice and back via periodic boundary conditions 



L{x) ^TrlP cxp 



I j Ao{x,xo)dx°j \ (6) 
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The static potential between sl qq pair is given by 
ln2?(a; — y) 



4. Duality @. 



Vqq{x-y) = 



(7) 



By cluster property, at large distances 
V{x-y) ~ C exp{~A\x - ylaNr) + K\{L)\'^ (8) 
A critical temperature Tc is observed such that 

forT<Tc {L)=0 V{r) = err (conf.) 

forT>Tc {L)^0 V{r) ^ (deconf.) 

r — ^oo 

Confinement is related to the presence of a linear 
potential at large distances: a is known as string 
tension. 

For SU{2) gauge theory the phase transition at 
Tc is second order, and belongs to the universality 
class of the 3d Ising model. Tc ~ 180 MeV§. 

For SU{3) gauge theory the transition is weak 
first order, and Tc ~ 270MeV§. 

In the presence of quarks the symmetry of 
which L is an order parameter, is not a symmetry 
any more. For massless quarks a chiral symmetry 
exists above some temperature Tc, which is spon- 
taneously broken for T < Tc, the pseudoscalar 
octet being the Goldstone particles. For mq ^ 
the chiral symmetry is again explicitely broken. 

It is qualitatively clear that confinement can 
produce a breaking of chiral symmetry: in a bag 
model chirality is inverted in the reflection on the 
confining wall. Numerical indications also exist 
that the two transitions take place at the same 
temperature. However the overall situation is not 
satisfactory. 

First of all if a phase exist in which color is con- 
fined, and a phase at higher temperature in which 
quarks and gluons are free particles, an exact or- 
der parameter should exist for this transition. 

In addition a strong theoretical hint exists that, 
when the number of colors Nc is sent large, with 
the constraint Ncg^ — A fixed, a limiting theory is 
defined, which does not differ much in its physical 
content from the realistic theory where Nc = 3. 

The expansion in 1 /Nc should be a convergent 
expansion [||. 

In this philosophy quark loops are non lead- 
ing 0{1/Nc) and therefore the physics of confine- 
ment, i.e. the symmetry, should be the same as 
in quenched approximation. 



Duality is a deep concept in field theory, in sta- 
tistical mechanics, in string theory. It applies 
to systems in d + 1 dimensions having non lo- 
cal topological excitations in d dimensions. Two 
complementary descriptions can be given of such 
systems. 

A direct description, in terms of local fields <&, 
in which topological excitations /i are non local. 
Symmetry is described by order parameters, the 
v.e.v. ($) of the fields $. This description is 
convenient in the ordered phase, or weak coupling 
regime. 

A dual description in which fj, are local opera- 
tors, and $ non local excitations. Symmetry is 
described by disorder parameters (/i). The dual 
coupling is gjj ~ This description maps the 
strong coupling regime of the direct theory into 
the weak coupling regime of the dual. Therefore 
it is convenient in the strong coupling regime. 

The prototype system[0 with duality is the 
1 + 1 dimensional Ising model: there the field $ 
is the spin variable a{i) = ±1. The dual con- 
figurations /i are 1 dimensional kinks. The dual 
description is again an Ising model in which the 
creation operator of a kink is = ±1 and the 
dual Boltzman factor is defined by the relation 



sinh2/3 



1 



sinh2/3D 



or /3 ~ 1//3_D. In the model T 
role of the coupling constant. 
Other examples are 



(9) 

1//3 plays the 



1) The 3d XY mode\^, which belongs to 
the class of universality of liquid He^. The 
topological excitations are abelian vortices. 

2) The Heisenberg magnet, where the topolog- 
ical excitations are 2d Weiss-domains [p3[ . 

3) the TV = 2 SUSY QCD|i| where the exci- 
tations are monopoles. 

4) Superstring Ad theories [p^. 

5) Compact U{1) gauge theories, where the 
topological excitations are monopoles [|6|, 



4 



The question is: what are the dual excitations 
in QCD? 

Two proposals exist in the literature, both due 
to G. t'Hooft|Jl|. 

a) Monopoles. Monopole condensation in 
the vacuum produces dual superconduc- 
tivity, and confinement of electric charges 
via Abrikosov flux lines (Meissner effect). 
Monopoles are defined by a procedure 
named "abelian projection" based on 
the choice of a local operator in the ad- 
joint representation, as discussed below. 
In a sense the mechanism is largely unde- 
fined, since there is a continuous infinity of 
choices for the operator used to define the 
monopoles and their interrelation is not un- 
derstood. A guess is that all these choices 
are physically equivalent iQ . 

The positive feature is that confinement is 
explained in a simple form: the chromo- 
electric field of a qq pair is channeled into 
an Abrikosov flux tube with energy propor- 
tional to the length, i.e. to the distance r 
between the pair jr^ , pOt 

V{R) = (jR 

a, the energy per unit length of the tube, is 
the string tension. 

2) Vortices. A vortex is a magnetic defect as- 
sociated to a closed line C . The operator 
B{C) which creates a vortex at some time 
t, obeys the following algebra, with the op- 
erator W{C') creating a Wilson loop along 
the line C 

B{C)W(C') ^ W{C')B{C)e^Y> (i^^27r) 

where ncc' is the winding number of the 
lines C, C . It can be shown that, whenever 
{W{C')) obeys the area law, {W{C')) - 
exp(— Areac), as the loop goes large, 
{B{C)) obeys the perimeter law, and vicev- 
ersa if {B{C)) obeys the area law, {W{C')) 
obeys the perimeter law. 

Area law implies that {L), the expectation 
value of the Polyakov line vanishes. One 



can define a dual Polyakov line L' wich is a 
vortex along a straight path crossing space 
from -oo to +00. Then for T < T^, (L) = 
0, for T > Tc, {L') = 0. {V} is a disorder 
parameter for confinement ||2l| . 

The definition of monopoles is associated to the 
choice of an operator ^(x) in the adjoint repre- 
sentation: we wfll speak of SU{2) gauge group for 
simplicity, but the procedure is easily generalized 
to SU{N). Let ^(x) = $(x)/\${x)\ be the direc- 
tion of $ in color space. A gauge invariant field 
strength Ff^,, can be defined as|Q 

f^,^^-g^,-^Hd^^ad,^) (10) 

= d^, - gAf,A 

The two terms in the definition are sepa- 
rately gauge invariant and color singlets: they are 
arranged in such a way that bilinear terms A^j^A^, 
or A^d,y^ cancel. Indeed 

F^, = ^d^A, - 9,4) - A dA) (11) 

The gauge transformation bringing $ in the same 
direction (say (0,0,1)) in color space is called 
abelian projection. After abelian projection 

F^, = d^Al - d^Al (12) 

The source of the dual tensor F*^ = l/2e^i,pcri^'"^ 
is a magnetic current 

Ju = d^P;, (13) 

and is conserved. The corresponding U{1) mag- 
netic symmetry can either be Wigner, and then a 
magnetic charge operator is defined and Hilbert 
space is superselected, or be Higgs broken, which 
implies the existence of at least one magnetically 
charged operator jj, such that {^) ^ 0. (/i) 7^ 
implies dual superconductivity. Notice that in a 
noncompact formulation j'^ = (Bianchi identi- 
ties). 

In principle the existence of dual superconduc- 
tivity can be investigated in any abelian projec- 
tion by looking at the v.e.v. (/i) of an operator 
creating a magnetic charge in that abelian pro- 
jection. 
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5. Construction of the disorder parameter 
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Once the dual topological excitations are iden- 
tified the disorder parameter can be constructed 
as the v.e.v. of their creation operator /i. The 
construction of /i[<i>] intcrms of the field of the di- 
rect description $ is an explicit realization of the 
duality transformation. 

The guiding idea goes back to ref. j27| and 
amounts to a traslation of the fields in the 
Schrodinger picture by the classical topological 
configuration. In the same way as 

e'P''\x) ^ \x + a) (14) 
= cxp ( z / SyT\{xo,y)(p{y- x) ] , (15) 



with n the conjugate momentum to the field $ 
[<^>{xQ,x),U{xQ,y)] = iS^{x~y), 
adds (p{x) to the field configuration 
/i(f)|$(x)) = |$(x) -f (^(f)) 

Adapting the above construction to a compact 
formulation of the theory as in QCD on the lattice 
is far from trivial, but has been donep^-p6t. 

A correct, gauge invariant, definition of the op- 
erator /Li for abelian projected monopoles ex- 
ists |^^. The resulting operator is a Dirac like 
magnetically charged and gauge invariant oper- 
ator, which can then have non vanishing v.e.v. 
without violating gauge invariance|p8[|. 

The construction has been checked on systems 
already studied and understood by other meth- 
ods, starting with the 2d Ising model|2^, where 
our operator fj, which creates kinks does indeed 
coincide with the dual variable cr* of ref. 1 1 1 1 . 

For the 3d XY model the phase transition was 
always conjectured to be produced by condensa- 
tion of vortices, but that condensation was de- 
tected in numerical simulations by observing an 
increase of the number of vortices in condensed 
phase: our construction shows that (/i) 0, or 
that the number of vortices is not defined in the 
vacuum of the condensed phase ]T^. 

For the Heisenberg model we have found [p^ 
that the Curie transition can be looked as an or- 



der disorder transition, and that the high tem- 
perature phase is ordered in the dual language, 
by condensation of non abelian 0(3) vortices. 

A numerical problem has also been solved; /x 
defined by eq. (|l^) is the exponential of an integral 
over the spatial volume: therefore it fluctuates 
typically as exp(L^''^), which is a wild fluctuation. 
In fact the disorder parameter (/i) is the ratio of 
two partition functions 



Z 
Z 



and has the fluctuations typical of a partition 
function. Z is the ordinary partition function of 
the theory, and Z — J corrisponds to a mod- 
ification of the action coming from the exponent 
in fi, eq.dl). 

A way out of this difficulty consists ||30| , p3|| in 
considering instead of (/i) the quantity 



(16) 



f3 



2N 



which contains all the relevant informa- 



tion of (/i). 

A typical shape of (/i) is shown in fig.l. The 
corresponding p is shown in fig. 2. The peak seats 
on the location of the phase transition. 
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Figure 1. (/i) for compact U{1). 
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6. Results for QCD||2|-||]. 
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Figure 2. p for compact U{1) 



The phase transition in principle takes place at 
infinite volume. The limit can be performed by 
a finite size scaling analysis as follows; if i5 is the 
critical index of fj, and r = {1 — T/Tc) the reduced 
temperature, 

then 



(17) 



The functional dipendence is dictated by the fact 
that $ is dimensionless: a is the lattice spacing, 
^ the correlation length, L the extension of the 
lattice. As r ^ 0, <^ goes large with a critical 
index i' 

a/S, can be approximated with (scaling limit), 
and L/^ can be traded with L^^^t. Hence 

and 



(18) 



pjL^^^ is a universal function of L^I'^t, i.e. it 
scales when plotted versus L^^^t. 

Since this scaling law is valid for the appropri- 
ate value of V and /3c, the determination of v, (3c 
and (5 is then possible. 



{jj) , or better p, has been studied for pure gauge 
SU{2) and SU{3) gauge theories across the de- 
confining phase transition, in a number of differ- 
ent abelian projections. 

A clear evidence has been obtained that, irre- 
spective of the choice of the abelian projection 

1) (/i) ^ in the confined phase (T < Tc). 

2) (/i) ~ exp(-fcLs), fc > for T > Ls be- 
ing the spatial extension of the lattice. This 
means (n) = at T > Tc in the thermody- 
namical limit. 



3) (m) 



1 - 



T 



4) Tc and v can be determined, together with S 
from the finite size scaling analysis sketched 
in sect. 3. 

Fig. 3 shows a typical form of p as a function of f3 
for different lattice sizes and SU{2) gauge theory. 

Fig. 4 shows the quality of the scaling equation 
(|r^), i.e. the validity of the extrapolation at infi- 
nite volume. 
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Figure 3. p for monopole condensation in SU{3). 
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Figure 4. Finite size scaling of p. 



Tc agrees with the determination made using 
the traditional order parameter (L), and so does 
V. For SU(2) V = 0.62(2) in agreement with the 
expectation that the transition belongs to the uni- 
versality class of the 3d Ising model, S = .20(8). 

For SU{3), V — 1/3, which means that the 
transition is first order, even if weak enough to 
allow a scaling region, and 8 = .50(3). 

The independence of the result from the choice 
of the abelian projection has been checked by 
performing the analysis in a number of different 
abelian projections, but also averaging on a very 
large number (infinite) of abelian projections]^. 

These results clearly indicate that confining 
vacuum is a dual superconductor. 

The same disorder parameter (/i) can be used 
in full QCD (including dynamical quarks) where 
it is well defined, contrary to (L) whose symmetry 
Zn is explicitely broken by the very presence of 
quarks and to (iptp), whose symmetry is broken 
by quark masses. 

The analysis for full QCD is in progress and 
of course requires a big computational effort. 
There the dependence on the masses of the quarks 
makes the finite size scaling analysis, eq.(18) more 
complicated. Preliminary results ||3l|] , however, 
indicate that also in full QCD there is a tran- 
sition from dual superconductor to normal across 
the transition. If these preliminary indications 



get confirmed by the quantitative analysis which 
is on the way, we would have a good order pa- 
rameter for confinement, even in the presence of 
(massive) quarks. In principle the two transi- 
tions, chiral and deconfinement, could take place 
at different temperatures; however all the existing 
indications are that they coincide. 

Moreover this would also reconcile confinement 
with the Nc — > oo limit. In that respect also an 
analysis of quenched gauge theories at A'^c > 3 is 
on the way, to perform the brute force check of 
the limit. 

Here again preliminary results show that the 
order parameter is weakly dependent on Nc at 
g'^Nc = const.^. 

However our results are still far from complete. 
They indicate that, whatever the dual fundamen- 
tal excitations are, they carry magnetic charge 
in all the abelian projections. A real theoretical 
breakthrough would be to identify such excita- 
tions, and to write an effective Lagrangian for 
them to describe QCD in the confined phase. 

Our results will hopefully help in solving this 
problem. 

I thank L. Del Debbio, M. D'EHa, B. Lucini, 
G. Paffuti, for discussions. Most of the work re- 
ported here is largely due to their collaboration. 
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